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Abstract 

We compute the exact partition function, the universal ground state degeneracy and 
boundary state of the 2-D Ising model with boundary magnetic field at off-critical tempera- 
tures. The model has a domain that exhibits states localized near the boundaries. We study 
this domain of boundary bound state and derive exact expressions for the function" and 
boundary state for all temperatures and boundary magnetic fields. In the massless limit we 
recover the boundary renormalization group flow between the conformally invariant free and 
fixed boundary conditions. 



1 Introduction 



Quantum field theory (QFT) in the presence of boundaries plays an important role in our 
understanding of a wide range of physical phenomena. Examples include various surface 
critical behaviors [0, impurity effects (e.g. the Kondo problem) p|, junctions in quantum 
wires dissipative quantum mechanics 0, etc. It is also important in open string theory 
1^ where the open string states are created by vertex operators built from boundary 
operators in the boundary conformal field theory. 

An important concept in boundary QFT is the boundary state. In Langrangian formu- 
lation the boundary condition may be expressed either in terms of a boundary action or as 
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a relevant perturbation of a conformal boundary condition. On the space-time of a cylinder 
of length L and circle length R with coordinate x running along the length and coordinate 
y running along the circle (Fig.l), the action may be written as: 

A= dxdyCbuik+ I dy^Bi{y)+ / dy ^bXv) (1) 

JV JBl JBr 

where $B; and are boundary fields on the left and right boundaries respectively. In 
Hamiltonian formulation there are two alternative approaches: (1) In the first one chooses x 
to be the space coordinate and y the Euclidean time coordinate. The Hilbert space of states 
Tif'^'' associated with a. y = constant time slice must satisfy boundary conditions Bi and 
Br at the left and right boundaries respectively. For example, in the L — > oo limit the states 
in Tif'^'' may be classified as asymptotic scattering states. Particles scatter with each other 
and off the boundaries. If the scattering is factorizable (admitting boundary Yang Baxter 
symmetry 0) then, relative to any boundary (and again in the limit of large L), one may 
choose either a basis of asymptotic "in" states or of asymptotic "out" states to span Tif'^''. 
The two bases are related by the factorizable scattering matrix. And the momenta of these 
asymptotic states are constrained by the boundary conditions in this Hamiltonian picture. 
The partition function is expressed as: 

Z = Trace^,B,Br(e'^^^) (2) 



'-L 



where is the Hamiltonian, and an n-point correlation function of local operators Oi(xi, yi), 
• • •, On{xn,yn) is givcu by: 

^ Trace^B.Br {Ty (Oi(xi, yi) . . . 0„(x„, yn) e"^^^)) (3) 

where Ty implies "y-ordering" . (2) In the other approach x is considered to be the Euclidean 
time coordinate and y the space coordinate. The Hilbert space I-Lr coincides with that of the 
theory with periodic (or antiperiodic) boundary conditions (i.e. the theory in finite space 
without boundary). The boundary conditions now appear as the initial state \Bi) at time 
a; = evolving to the final state \Bj) at time x = L. These states are called boundary 
states and encode all information about the respective boundary conditions. The partition 
function is now expressed as: 

Z = (5,|e-^^«|5,) (4) 



where Hr is the Hamiltonian in this formulation. An n-point function is expressed as: 

{Br\%Ol{Xi,yi) . . .On{Xn,yn)\Bl) 
{Br\Bi) 



(5) 



where 7^ now implies "x-ordering" . 

The two dimensional Ising model exhibits a rich variety of phenomena many of whose 
characteristics can be calculated exactly. It has often formed the basis of our theoretical 
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understanding of phase transitions and critical phenomena, the renormahzation group and 
general quantum field theory. The model on a manifold with boundaries and with a boundary 
magnetic field has several interesting features. The boundary model is integrable and its 
exact boundary scattering matrix is known for all temperatures 0. At critical temperature, 
i.e. in the massless case, the model describes a renormahzation group (RG) fiow from the free 
to the fixed conformal boundary conditions [|1^, and exhibits a monotonically decreasing 
universal function" 0] akin to Zamolodchikov's c function [|12| in the bulk. This g function 
and the full boundary state were computed exactly in ||13[. They have also been computed in 



14| using the thermodynamic Bethe ansatz (TBA) approach for boundary theories [|l5l , |T 



In general, the boundary TBA approach yields only the ratios of the g factors due to some 
difficulties mentioned in |[T^ . In this paper, we compute the exact partition function, ground 



state degeneracy factor g and boundary state for the massive case using the approach in 13 



2 2-D Ising Model with Boundary Magnetic Field 

The 2-D Ising model with boundary magnetic field has been a subject of investigation for 



quite some time. In [0, (see also [|^) the model was studied on the lattice and the 
boundary contribution to the free energy and the boundary spin spin correlation function 
were computed (see also pO|). The free (i.e. boundary spins are free) and fixed (i.e. all 
boundary spins are fixed to +1 or —1) limits of the boundary condition were studied in the 
context of boundary conformal Field theory and the boundary spin operator was identified 
in . At critical temperature, the boundary RG fiow from the free down to the fixed 



conformal boundary conditions was identified in 0. The model was studied in the context 
of boundary integrable field theory in 0. In ||21| the local magnetization was computed 
exactly in the massless theory using operator product expansions and were found to be 
hypergeometric functions. Differential equations for the local energy and magnetization 
have been derived recently in The exact partition function and boundary state were 
computed for the massless theory in [p!3| , 

As is known, the continuum theory of the Ising model is described by free Majorana 
fermions. On a manifold T) with boundaries Ei, i = 1 . . .n, the model is described by the 
action: 

A = — / dx dy [ip dz ip + ip dz ip + imipilj] 
27r Jt> 

+ E{ / dt[-^^^ + ^ad]+ihJ dta{t){e]^ + e]tlj){t)\ (6) 
j=i \ ^ J / 

Here z = x + iy and z = x — iy are the usual complex coordinates on the manifold V; 
z = Zj{t) and z = Zj{t), with t being a real parameter, parametrize each boundary Bj] 

ej{t) = ej{t) = ^ are components of the tangent vector {ej , ej) to the jth boundary 
with ej{t)ej{t) = 1; ip{x,y), il){x,y) are the free Majorana fermion fields and m is the mass; 
a{t) = is a fermionic boundary field (see P, |2T|), anticommuting with ip and ip, with the two 
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point function: 

{a{t) a{t'))free = lstgn{t-t') (7) 



2 

and a = ^ ; h is the appropriately rescaled external boundary magnetic field with the 



dimension of [length] 2 . The integrand in the last term in (P) is the boundary spin operator 
as (see 1,0, ig): 

aB{t) = a{t){e]tlj + e]4j){t) (8) 

The first two terms in comprise the bulk action and the fermion mass m is related to the 
temperature by 

m ~ (T, - T) (9) 

As is known, the field theories corresponding to the ordered (T < Tc) and disordered (T > Tc) 
phases are related by the duality transformation ip ^ ip , ip ^ —ip and are equivalent. We 
will assume the field theory for the low temperature phase i.e. m > throughout this paper. 
From ip, ip satisfy: 

d d 

{— + i—)'^ + imtp = (10) 
ox oy 

d d - 

^d^~%^'^~''^'^^^ ^^^^ 



in the bulk and 

at each boundary Bj, where 
and 



;A+,A)^^.(t) = (^_,A)V;,(t) (12) 



A = 47rr (13) 



V^,(t) = e]it)^{Z^it)) , ^P^it) = e|(t)^(Z,(t)) (14) 
For our cylindrical geometry of Fig.l boundary conditions (|12]) take the form: 

(- — h iX)uj'ip{0, —y) = (- iX)uj'ip{0, —y) at the left boundary a; = 0, and (15) 

cLy ay 

(- — h iX)uj^jj(L,y) = (- iX)iu^jj(L,y) at the right boundary x = L (16) 

dy dy 

where u = and uj = e~ . 

In the Hamiltonian picture (which from now on we call the "L channel") where x plays 
the role of Euclidean time and y is the spatial coordinate, ip, ip admit the following decom- 
positions in terms of plane waves in the Neveu Schwarz (NS) (where ip, ip are antiperiodic 
along y) and Ramond (R) (where ip, ip are periodic along y) sectors: 

'iP{x,y) = J2 [a{9n)u;e'-^e-'^^''°'^^-''+'''''^'^-y'^ +a\9n)ooe^e'^'^'''''^^^^^^ (17) 
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where 



and 



27r 1 

sinh6'„ = — —{n H — ) , -n, G Z in the NS sector 
mR 2 



27r 

sinh 9n = — —n , n G Z in the R sector 
mR 



(18) 

(19) 
(20) 



In (|T3), (IH) o-iOn) and a^(6'„) are the annihilation and creation operators for a Fermi particle 
moving along the circle with momentum m sinh 6^ and energy m cosh 6^, On being its rapidity, 
and satisfying standard anticommutation relations. A boundary state \Bh) at a; = satisfies 



^ a \ x,^n (Q - cosh e, 

a(-^„)-.tanh- (^^^ 



\Bh) = 



where 



Q 



A 



-11 



v. 



(21) 



(22) 



m m 

and a similar equation holds for a boundary state at x = L. These are obtained by imposing 
(|15D , (|16D on {\[7\) and ([T8|) . Notice that the factor in the second term in ( pi]) is precisely 
the boundary scattering matrix in the cross channel /C(6'„) = 7^(y — 9n) discussed in 0. 

In the cross channel (henceforth we call it the "R channel" ) in which y is the Euclidean 
time and x is the space coordinate, the mode expansions for ip, ip are: 



— ^ ^—m(i sinh Sjx+cosh diy) 



^{x,y) = J2 



—i b{6i)e 2 e 



-i- iri(i sinh 9ix+cosh Siy) 



p—m^i sinh Sjx+cosh ^jy) 



«)e2 e 



(23) 



(24) 



for y G [0, i?) and 

ip{x,R) = —ip{x,0) and ijj{x,R) = — ^/'(x,0) in the NS sector, (25) 

'?/'(x, R) = ip{x, 0) and i'ix, R) = i^i^x, 0) in the R sector (26) 

Here b{6i), h\9i) destroy and create fermions (in the R channel) which now move along x 
with rapidity 9i. These fermions scatter off the boundaries with the scattering matrix TZ{9) 
given by [Q: 

h\-e) =n{e)h\e)- (27) 

9 Q + i sinh6' 



n{e) = i tanh 



(28) 



4 2' Q -2 sinh ^ 

which can be obtained from (plSj), ( p3D and (p^ . The rapidities 9i in ([23|), (p^) are solutions 

2 



of: 



1 + X = 



where X 



,,2imL sinh S 



,tTC 



tanh (— 



^, Q-isinhO 



4 2'\g + isinh^ 



(29) 
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which we get from ([TB|), (PB]), and We notice that although 6^^ = is a solution 



of (p9|) it is unphysical since 7^(0) = —1 implying that the wave function for this zero 
mode vanishes. In the next section we start with the partition function (|]) as a trace over 
asymptotic states (p3D, (0), p9| ) and make explicit transformations so as to express it as 
a sum over states in the cross channel. In the large L limit this will yield us the boundary 
state. 

3 Partition Function 

As mentioned in the earlier section, the spinor fields if), '4) can be either periodic or antiperi- 
odic along the circle of the cylinder. We have been calling these the Ramond (R) and the 
Neveu Schwarz (NS) sectors respectively. In the R channel the states associated with the 
?/ = and y = R lines are identical in the Ramond case while for the Neveu Schwarz case 
they differ by a phase of vr. The partition function (0) in this channel generalizes to Z+ and 
Z_ in the Neveu Schwarz and Ramond sectors respectively, where: 

Z± = Trace^B.B^ ((±1)^6"^^^) , (30) 

T being the fermion number operator. For our boundary Ising model (P) the partition 
function ( pOf ) is given by: 

Z± = e-^^^ n (1 ± e"'"^™^^^') (31) 

where Ej^ is the ground state energy and Qi 's are solutions of (|29|) . The free energy F-j. is 
given by: 

- RF^ = -REl + log(l ± e-"*^™^^^') (32) 



which can be written as (as in [|T3]): 

- RF± = -REl + — I d9 log(l ± e'^^™^^^) (33) 

2711 Jc 1 + X 

where X' = ^ and the contour C in the complex 9 plane is shown in Fig. 2 corresponding 
to the case when Q > 1. We note here that X{-e) = 1/X{e) and X(z7r - 6) = X{0) The 
positions of zeroes and poles of 1 + X in the complex 6 plane fall into three distinct domains: 

(i) Q > 1 : In this case all zeroes of 1 + X lie on the real axis and they cluster closer together 
further from the origin and with increasing L. There are no poles in the upper half plane in 
the physical strip. The picture is shown in Fig. 2. 

(ii) < Q < 1 : In this domain also all zeroes are on the real axis as in (i). There is a pole 
on the negative imaginary axis at —i sin~^ Q in addition to the ubiquitous pole at —i'n/2 but 
none in the physical strip of the upper half plane. The analytic structure is shown in Fig.3. 

(iii) Q < : This domain is perhaps more interesting. There is a pole on the positive 
imaginary axis at 6' = — i sin^^ Q and possible zeroes on the imaginary axis as well. The pole 
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corresponds to a boundary bound state in the R channel. We will discuss this domain in 
greater detail later in the section. 

Let us first consider domains (i) and (ii) where all solutions of (|29|) are real and 1 + X has 
no pole in the physical strip. 



3.1 Q > : Domain of No Boundary Bound State 

Here we treat domains (i) and (ii) together. First we compute El : 

1 If X' 

Er = — m cosh 6'/ = / d9 m cosh 6' (34) 

2^^o AmJc 1+X ^ ' 



We can rewrite the integral on the right above as: 



El = — —if de + f d9 - I rf^l^^mcosh^ 



(35) 



where the contours C+, C_ and Co are shown in Fig. 4 and Fig.5. The last integral in (^) is: 

/ de -cosh^ = — 36 

SiriJc. l + X 4 ^ ^ 

We rewrite the C_ integral in (BH) as: 



— / de cosh e = / de — - — ^ h / de — cosh e (37) 

MiJc. l + X SniJc- X^{1 + ^) 8m Jc- X ^ ' 

After changing variable ^ — — ^ in the first term above and expanding the second term we 
obtain: 

El = / de cosh^ / ci^cosh^e / de\l + —^ ^ (38) 

4 A-niJc+ l+X 2t^ Jo 27r io ^ Q^ + sinh^^^ ^ ' 



The third term in (p8D contains the nonuniversal bulk intensive free energy which we set to 
zero. We recognize the fourth term to be the boundary intensive free energy and denote it 
by e\ : 

m ,^ r, 2Q cosh^ e , 

We now consider the second term in (pH]). We shift the contour up by y (see Fig. 6) (note 
that we can do this trivially since the integrand has no singularity in this region), change 
variable ^ — (y — 0) and integrate by parts to obtain: 

I f X' I f°° 

— de -cosh^ = — / c/^cosh^log(l + yL(^)) (40) 

AmJc+ l + X An J-oo 
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where 



-2mLcoshe ^^^-^2 ^ I Q COsh 



2 VQ + cosh^ 



Finally we combine (|8D, (H), (|0|) to obtain: 

171 171 

EL = - + 2e^-— d^cosh^log(l + YL{e)) 
4 47r j-oo 



(41) 



(42) 



Now let us compute the second term in (^). As we did for E^, we rewrite the integral 



as: 



RF± = -REl + — / de 
Am \Jc+ 



where 



YR = e 



de 



-mi? cosh 9 



d9 



X' 



Co 



1 + X 



log(l ± Yr) 



The last integral is trivial: 



de-^ioR(i±YR) 



lofffl ± e 



-mR\ 



ATxiJco l + X"°'~ 2 
Once again, in the integral along contour C_, we rewrite in the integrand: 



(43) 
(44) 

(45) 



X' 



X' 



X' 



1 + X X2(l + ^) + X 



(46) 



Then we transform variable ^ —0 as we did in (p8| ) and get: 



— / de 



X' 



Am Jc- l + X 



logfliFR 



/ de 



X' 



+— / de mL cosh e log(l ± Yr) + — / 

27r j-oo 27r J-i 



de 



Am Jc+ l + X 
1 ^ 2Qcoshe 



logfliF, 



R) 



cosh e + sinh^ e 



log{l ±Yr) (47) 



We recognize the second term in (^) to be the Casimir energy: 

n(±) 



1 roo 

E'^' = / rfemcosh^log(l ± e-™^^°^^^) (48) 

27r J-oo 



where the ± signs, as usual, correspond to the NS and R sectors. Next, for the first term in 
(^81) , we apply a Wick rotation of the momentum axis in the complex plane by shifting the 
contour C+ in the e plane up by i7r/2 (see Fig. 6). This enables us to go to the cross channel. 
Then we change variable 6* — > (6^ — y) and get: 



— de- 



X' 



2m Jc+ l + X 



logfliFfi 



1 



2TTi Jc 



de- 



Y[ 



l+Fr 



logfl ± X 



R) 



(49) 
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where the contour C' is shown in Fig. 7, Yl is given by ( |1T| ) and 



Xnie) = Yn{- -e)=e 



imRsinhO 



After integrating by parts (^) yields: 



de 



X' 



log(l±F«) 



1 



X' 



R 



■iog(i + ri) 



27ri Jc+ l + X °' ^ 27ri Jc' Xr ± 1 
We rewrite the integral along C as a sum of integrals along and (see Fig. 7) 



(50) 



(51) 



— I de ^^^^ iog(i + Yl) = — { I de + I 

Jc' Xr±1 ' 2Tii \Jc'^ Jc'_ 



2m Jc' X 



de 



Xr±1 



log(l + Yl) (52) 



Our aim is to express the free energy as a sum over states in the cross (i.e. L) channel. So 
with hindsight we change variables e —e in the C_ integral in (|52|) and obtain: 



— / de 



X' 



R 



2m Jc' Xr ± 



ci^cosh^log(l + yL) 



mR 

Jo 



Thus (Ell) becomes: 



de ^^\og{i± Yr) 



2Tii Jc+ l + X 



de- 



X' 



R 



2m Jc" Xr ± 1 
mR 1-°° 



iog(i + ri) 



271 



POO 

/ rf^cosh^log(l + 7^) 

JO 



where C = C.-\ — C_ is shown in Fig. 8. The first term above evaluates to: 



de 



X' 



2m Jc" Xr ± 1 



log(l + Yl) 



2n 1 

logS+ = J2 log(l + MuJi)) , msinhcjz = ^ ^ Z + - for (+) 

i>o ^ ^ 

27r 

logS_ = ^ log(l + iL('^n)) , msinhciJn = —n , n G Z for (— ) 

n>0 ^ 

Finally combining (||), (||), (||), (0), (|8|), (H) and (||) we obtain: 



- = ^ - 2£Ai? + log S± - - log(e— ± e"— ) 



1 /-oo 

27r j-oo 



1 ^ 2Qcosh^ 



cosh^ g2 + sinh^^ 



log(l ± e" 



-mi? cosh Q 



(53) 



(54) 



(55) 



(56) 
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We identify the last two terms as the universal ground state degeneracy: 



1 /""^ 
log g±{R) = — de 
47r J-oo 



1 ^ 2Q cosh 9 



cosh^ Q2 + sinh^^ 



logfl ± e" 



-mR cosh d\ 



--log(e— ±e- — ) (57) 

Thus from (^) and (|57D the partition function is: 

Z± = e-^^«*'-2-^^S±(?2(^) (58) 

Expressions ( ^6|) and ( ^8l) with ( ^5|) and (^) are exact. Notice that in the L ^ oo limit 
Let us now look closely at the Q < domain. 

3.2 Q < : Domain Containing Boundary Bound State 

In the Q < domain, i.e. for sufficiently weak boundary magnetic field, the boundary 
scattering matrix has a pole in the physical strip at ^ = —ism~^Q. Correspondingly 
1 + X has a double pole at this same 6. A careful analysis shows that, for sufficiently large 
distance L between the boundaries, there are also two zeroes pinching the pole from above 
and below and approaching it along the imaginary axis exponentially fast with increasing 
L. There are also two more zeroes on the negative imaginary axis located symmetrically to 
the upper ones and approaching 9 = isin~^ Q at the same rate. There is however no pole 
there. These zeroes correspond to states (in the R channel) localized near the boundaries. 
In fact there are two such states because the particle can be bound near the left or the 
right boundary and the zero above the pole corresponds to a symmetric (S) superposition 
of states localized near the left and right boundaries while the zero below corresponds to 
their antisymmetric (A) superposition. Each boundary may be looked upon as an infinitely 
heavy particle sitting there with the light fermions scattering off it. At weak boundary 
magnetic fields they form bound states - the fermion gets trapped by the heavy boundary 
particle. In the limit of large L the S and A superposed states are degenerate (each bound 
fermion does not see the opposite boundary and their wave functions do not overlap) while, 
for finite L, they split by an amount exponentially small in L due to the overlap between the 
left and right localized states (the bound fermion can now tunnel through the finite barrier 
and become bound to the other boundary particle). Moreover, when L decreases one zero 
each from above and below the real axis, which are closer to the origin, eventually leave the 
imaginary axis and go to the real axis joining their friends there. The picture for sufficiently 
large L is shown in Fig. 8. As we will see our results will not depend on L. 

So now we include the contributions from the extra solutions to (E3) in the integrals in 



(H and (H). We have: 

X' 



Er 



Airi 



f de + f de + f 



de 

0"1 ■Jo-2 



l + X 



mcosh^ (59) 
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and 



-REt 



1 

27Ti 



dd + dO + 



de 



X' 



1 + X 



log(l ± e 



-mR cosh d 



) (60) 



where ai and (T2 are small closed contours enclosing the zeroes below and above the pole 
respectively in the upper half plane. We repeat the same steps as before in computing the 
integrals along contour C. However, when we shift the integration contour C+ up by y in 
obtaining (^OD and p9|) we have to avoid the singularities of in the integrand at the two 
zeroes and the (double) pole of (1 + X) (see Fig. 9). In the process we pick up contributions 
from each of these singularities (see Fig. 10): 



, dd- 

2'7ri Jc+ 1 



X' 



X 



(...) 



1 

2Td 



de 



de 



de - 



0-2 



de 



pole 



X' 



X 



(■■■) (61) 



Here C_|_ is the same as in ( pSf ) and (^). The last two integrals cancel exactly with corre- 
sponding integrals in (^) and (|60|) and we are left with the integral along apoie only. Thus 
the extra contribution, in effect, comes only from the pole, irrespective of the number of 
zeroes on the imaginary axis (of course we have to include all these zeroes in (|59D , to start 
with). In the light of the above discussion, our result is therefore independent of L as long 
as we include a a contour for each zero (of 1 + X) in the physical strip. Thus we have: 



1 



pole L -r yy 



m cosh e 



-mR cos u 



where 



u = — sm 



■'Q, g<o 



(62) 



(63) 



Expression (^) gives the bound state contribution to the boundary energy. For Q < the 
intensive boundary energy (let us now denote it by e\) is therefore: 



24 = -f rde[i + 

ZTT JO 



2Qcosh^^ . 
Q2 + sinh2^ 



— mR cos u 



And we have: 



1 

2T(i 



de^ 

pole 1 + X 



log(l ± e" 



-mR cosh I 



21og(l ±e- 



Q <0 



-mR cos u\ 



(64) 



(65) 



implying that now in the presence of boundary bound state in the R channel, the L channel 
ground state degeneracy function (which we denote by g\. for Q < 0) is: 



9l{R) 



1 

47r 



de 



1 



2Q cosh^ 



log(l ± e 



-m/?cosh 9 



■ log(e 



cosh^ Q^ + sinh^e_ 

) +log(l±e-™^^°'") 



mR mR 

2 ± e 2 



g <o 



(66) 
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In an alternative interpretation, we can consider g'^ as an analytic continuation of g± given 
in (|57|) with the integral in ( ^Tf ) now being along the contour C^^^ shown in Fig. 11 instead 
of along the real axis. The contour surrounds the two poles of the integrand at ±i sin~^ Q. 
Notice that the term in [■ ■ ■] in the integrand is precisely the derivative of the phase shift 
from boundary scattering and appears also in the boundary thermodynamic Bethe ansatz 
approach [|16|. g± defined this way covers the whole domain of Q and we treat this as our 



general definition for g±. Denoting this by g^\ we have 



^ ' 47r7c{s) 



2Q cosh e 



cosh^ Q2 + sinh^^ 



log(l ±e' 



-mi?, cosh S\ 



1 , raH. mH. , , , 

- log(e 2 ±e 2 ) (67) 




g > 
g <o 



(68) 



where the contour C^®^ is shown in Fig. 11. In the same way, we define the intensive boundary 
energy e^^^ for all Q: 



2e 



(Q) 



-^f de[i+ y^^^h^f ] 

47ric(s) ^ g^ + sinh^^^ 



(69) 



AQ) 



g > 
, g<o 

The partition function in the whole domain of Q is therefore given by: 



(70) 



(71) 



We will now compute the boundary state in the next section. 



4 Boundary State 

We consider the low temperature phase T < of the Ising system as we have been doing. As 
we discussed earlier the boundary state belongs to the L channel Hilbert space Hr. In this 
space there are two vacuum states |0) and \a) lying in the Neveu-Schwarz and Ramond sectors 
respectively. In the infinite volume (large R) limit the allowed momenta are continuous. 
The ground states |0) and \a) become degenerate. They are respectively the symmetric and 
antisymmetric combinations of the vacuum states |0+) and |0_) which are the broken Z2 
symmetry ground states and correspond to all spins pointing up and all spins pointing down 
respectively: 

|0) = 4s ( |0+) + |0_) ) , W) = ^{ |0+) - |0_) ) (72) 
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If n is the overall spin flip operator, we have n |0) = |0) and 11 |(t) = 
and Z_ represent contributions from the 11 = 1 and 11 = —1 respectively. We expect that: 

Zhh ± Zh^h 



where 

implying that 



^hh 



(5^|e-^^«|fi^) and = e'^^^ 



Zhh 

For large L we have: 



and Zh-h = -{zf^ 



(+) 



{-) 



^hh 



'h-h 



-LE 



(+) 



mBh)?-e 



-LE 



{-) 



{{o\Bh)f 



a). Obviously Z+ 
at: 

(73) 
(74) 
(75) 

(76) 
(77) 



where we have used {0\B_h) = {0\Bh) and {a\B_h) = and we have chosen the phase 

of |0) to be such that {0\Bh) and (cr|-B/i) are real (see |^). Thus, from ([75|) - ([77|) it follows 



that: 



(0|5,,) = ^e-«^A-^^:^^(i?) and {a\Bh) = 



(Q) 



,(Q) 



-Re 



(Q) 



V2 



V2 



(78) 



As was shown in |^, the presence of an infinite set of integrals of motion for a system with 
boundary makes the boundary state take a simple form with contributions from pairs of 
particles of opposite momenta (and possibly a zero mode) only. Following |^ and as was 
done in [|l^ for the critical case, we write down the boundary state from (plTf ), (|l^) and (iTSf ): 



-Re 



(Q) 



gf\R) exp \ Y: \ |0) 



.n=0 



± g^3\R) exp \ JCiOn) a^e^) at(-^^„) \ \a) 



(79) 



where sinh6'j = g^\R) and e^^^ are given by ( |67|) and (|69| ) respectively, and from 



/c(^) = 7^( 



ITC 



i tanh 



^ / Q-coshg ^ 
2 I g + cosh0 



(80) 



5 Asymptotic Limits 

We can take the (bulk) ultraviolet limit mR — in two ways: (i) we can either make m ^ 
i.e. increase the temperature T ^ on the same cylinder, or (ii) we can reduce the size 
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R of the system R ^ while maintaining the same temperature so that the system size 
becomes much smaller compared to its correlation length. 

Let us discuss case (i) first. As m — > 0, Q increases, the boundary bound state, if initially 
present, disappears and = g±. The theory now describes a renormalization group flow 
in the space of boundary interactions interpolating between the free and the fixed conformal 
boundary conditions 0, The massless particles group into right and left movers which 
do not see each other. We have to obtain the mR limits of 71(0) and }C{6) carefully. 
We do this by shifting the rapidity occurring in TZ{6) by infinite amounts: 6 — > log ^tb + ^ 



for the right movers and 9 



mR 

log — 9 for the left movers. For the continuation of TZ 
in the cross channel, i.e. /C(^), however, we must shift 9 — >■ — log^ — 9 for right movers 

IC{ik), k 



and 9 log:^ + 9 for the left movers instead to obtain the correct /C(/c) 
being the momentum of the massless fermions (e.g. for right movers k 



1 

R' 



Moreover, 



in the mR 



limit the integral in (^) falls into three regions: (i) 9o ^ log^, i.e. 
h'^R < 1 where 9o = sinh"^ y/Q'' - 2 is an extremum of $(^) = , (n) 9o - Tog; ^ 



mi?' 



i.e. h^R ~ 1 and (m) 9o 
(ii) the following hold: 



^ log-^, i.e. h'^R 3> 1 It is easy to see from (|57|) that in region 



9+iR) 



r(« + i) 



21/4 



a 



with regions 
tively: 



and 



9-iR) = 
where 

above appearing as h^R 

- = \p2 and 
and 



a 



r(a + 1) ve. 
2h^R 

and h^R 



oo limits of 



9- 



1 



9- 





1 

2^3 



for h^R- 
for h^R 



0, 



oo 



(81) 

(82) 

(83) 
respec- 

(84) 
(85) 



We then obtain the boundary state in the massless limit from ([7UD, ( P^ and 



-Rex 



a 



1 



T(a + 



■J- exp < I 



± 



21/4 



E 



n 



2 - a + 
^ -a' 



r(a + 1 



a \ . y-^ n — a 
expii 2^ 



T{n + a 



|0) 



1^) 



^6) 



sions (B2D and 



where creation operators for right and left movers of momentum kj = ^j. Expres- 



were obtained directly in [|T^, |1J] in the case of the critical Ising model. 
As was shown in [jl3|, at the conformal limits, ( 186D yields for: 

oo 



\Bfree) = exp(i 



oo 

71+ i r. 

n=l 



.i)(|0)-k)) 



7) 
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where |e) is the energy density operator in the bulk conformal field theory and we used: 

at a\ |0) = i\e) ( 



and for 
h^R^ oo 



oo 

|5/ixed±) = -^exp(-i l]ai+ial+i)(|0) + \e)) 



n=l 



oo 



^ ^ exp(-i ^ (89) 

^ n=l 

Both ( |87| ) and ( p9D agree with Cardy's results ||10| 

Now we consider the case (ii) way of taking the mR — > limit, namely we let i? — > 
while keeping the temperature constant. For finite Q it is easy to see that as i? ^ the 
system reaches the conformally invariant limit of free boundary conditions. We readily see 
that in this limit ( |67D yields g+ = and = for both Q > and Q < domains (here 
we note that for Q < the second term in the integral in (^) yields log instead of log \/2 
that one obtains for Q > 0) and expression (^) for the boundary state. 

niR oo 

Now we consider the case mR oo when the size of the system is much larger compared to 
its correlation length. From ( pTj) we readily see that in this limit, 

9'^\R) = 9^^\R) - e-^ (90) 
From ( p7D and ( [7D| ) it is easy to obtain : 

\B±h) = ^=(|0) ± |cr)) + excited states (91) 
v2 

as we expect since these are just the broken symmetry vacua |0+) and |0_) (^) discussed in 
the previous section. 



6 Conclusion 



Expressions (67) for the "g function" and (|79D for the boundary state are the most important 
results of this paper. We have computed these quantities by starting with the partition 
function expressed as a sum over states (satisfying the boundary conditions) in one channel 
and making explicit transformations so as to express it as a sum over states in the cross 
channel. The domain Q < 0, in which there are states localized near the boundaries (in the 
R channel), is particularly interesting. We found single expressions for the g function and 
the boundary state valid in all domains of Q by analytic continuation, whereby the contour 



15 



of integration is made to surround some relevant poles. It would be interesting to see if 
this could form a general prescription for treating boundary bound states in obtaining the g 
function in general integrable quantum field theory with boundary, and it would be desirable 
to gain a proper physical understanding underlying this treatment. 
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